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We determine doubly heavy baryonic potentials as a function of the distance between the two 
static sources, coupled to a light relativistic quark, for different quantum numbers. We use the 
variational method to compute the ground state and the first two excitations. These can be used 
as an input to nonrelativistic models or to NRQCD calculations of properties of doubly heavy 
baryons. We compare our findings with a factorization model. We employ all-to-all propagator 
methods, improved by an additional hopping parameter expansion and Wuppertal smearing on 
N f = 2 QCDSF configurations. 
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Figure 1: The different models for the static-static-light baryon. 



1. Introduction 



The coupling of quarks inside a doubly heavy baryon is debated and there are various models 
on the market like A- or 7-fiux tubes between all three quarks or the formation of a diquark which 
couples to the third quark. The ground state doubly heavy baryon potential has been studied pre- 
viously on the lattice [1], Ref. [2] suggests that a diquark formation of the two heavy quarks is 
dynamically favored and that the light quark perceives this as an almost pointlike object. Therefore 
a QQq baryon with two heavy quarks is related to a heavy-light Qq meson [3]. In the static limit 
where the spin of the heavy quark completely decouples, the main difference between these two 
systems is that the QQ diquark is spatially extended. In this work we check up to what QQ separa- 
tions this scenario holds. In figure 1 we depict different situations. We ascribe a phenomenological 
Compton wavelength A ~ \/m q to a light quark of mass m q . For r <C A~ ! the light quark cannot 
resolve the two static quarks (left picture). When either the light quark mass or the static quark 
distance r is increased, the light quark can discriminate between the static color sources and can 
"choose" its preferred localization. Symmetry considerations suggest the scenarios where the light 
quark sits on top of one of the static quarks or in the middle between them to be of particular inter- 
est. The first case corresponds to Qq diquark formation while in the second case no diquarks are 
formed. Unlike in the physical situation of finite heavy quark masses where the average distance 
(r) is dynamically determined, in the static limit we can change r to interpolate between the QQ- 
diquark picture at small r and the other more involved scenarios. At r = the spectra of the QQq 
and Qq situations are (numerically and analytically) identical. As r is increased the energies will 
increase. For instance, Vqq = Vqq/2, up to a constant, at least to the first two orders of perturbation 
theory and in the string picture. We denote mass levels of the mesonic Qq system by m-Q q . Thus, in 
the <2<2-diquark picture, that at least should hold at r < A 1 , we expect the factorization, 



V QQq (r) 

which can graphically be depicted as, 
exp 



1 

m Q q + 2 V QQ 



(1.1) 




exp 




(1.2) 



where straight lines denote gauge transporters (including the static propagators) and wiggly lines 
stand for light quark propagators. The QQq potentials only depend on the distance and a D'^ 
representation (see below). The light quark and gluon degrees of freedom have been integrated out. 
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One can use the resulting potentials as an input to a Schrodinger equation, in the heavy quark limit, 
and apply the formalism of NRQCD for corrections. 

2. Simulation details 

We use 68 QCDSF configurations [4]. The specifications are displayed in table 1. We use 



Lattice points 




^sea 


P 


Csw 


m K [GeV] 


a [fm] 


L [fm] 


16 3 x 32 


0.1355 


0.1355 


5.29 


1.9192 


0.783(11) 


0.084(1) 


1.3 



Table 1: Wilson-Clover N f = 2 lattices. 



all-to-all propagators with 300 stochastic estimates on each configuration. We apply the hopping 
parameter expansion [5, 6] to reduce the stochastic noise. APE-smeared [7, 8] gauge links were 
used for the spatial gauge connectors in the Wilsonloops and the link covariant displacements. For 
the static quark propagator we use temporal links that have been stout smeared [9] once, to reduce 
the static energy and thus to improve the signal to noise ratio. The quark fields of the all-to-all 
propagators are Wuppertal smeared [10] using APE-smeared gauge transporters and we applied 
the variational method [11, 12] to extract the groundstate (GS), the first excitation and second 
excitations (FE, SE). With this setting we follow the methods described in [6]. We compute the 
baryon correlator with two static quarks in various geometries. The light quark can be in the middle 
which we will refer to as □, 

I ? I, (2.1) 

or the light quark can be placed at one of the static quark positions, symmetric (©) or anti- 
symmetric (0), 

(Q-OMH-0) « 

All these correlators are realized by using the appropriate gauge connectors & in the baryon inter- 
polator 

0(x,x',x") a = e abc Pf Q{x)% 0>{x,x") a , a (Q{x%&{x\x") m T^q{x")l) , (2.3) 

where all coordinates x are at the same time t , P + = I (1 + 74) is the projector to positive parity and 
r is one of the operators from table 2. We denote the separation of the two heavy quarks by r. We 
use the same operator T for source and sink so that the correlator we compute is characterized by 
r,r, the temporal extent of the correlator T and the geometry ( □ ,©,© ). For r = 0, CD and © are 
degenerate. 

The octahedral group Oh is the cubic group of rotations on the lattice, with the addition of 
parity. Its irreducible representations describe point-particles at rest on a discrete lattice. The map- 
ping between these and the continuum quantum numbers is not unique and so an lattice irreducible 
representation can correspond to multiple continuum / states. 
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In our case, at r > 0, this symmetry is broken to the cylindrical subgroup C 0] x . In the contin- 
uum this corresponds to D^, C 0(3). As we are interested in baryons we have to take a spinorial 
irreducible representation (irrep). We refer to these representations with half integer values of J as 
0' h or D'^ h ID' Ah . In table 2 we list the operators that we have used to calculate different correlators 







r = 


r >0 


Wave 


Operator 


0'(3),0' h 




S 


75 


1+ r+ 
2 ' °i 


2g> G ls 


P 


1 






P+ 


7iA] - 72 A 2 cyclic 


\ ' H 


2» II' G 2» 

L - 1 ' Gi » 


D 


75(71^1-72^2) ©cyclic 


3+ H + 
2 ' 


l g II' G 2g 
5g - 1 ' G l§ 


D + 


71 A 2 A 3 + 7 2 A 3 A! + 73A1A2 


5+ r + 
2 ' G l 


2g' °ig 


F 


75(71 A 2 A 3 + 7 2 A 3 A! + 7 3 AiA 2 ) 


1 > G r 


iAu> Gi « 



Table 2: List of operators and representations. 



and their corresponding quantum numbers. We denote the link covariant displacement as, 

A^ri(x) = U^(x)ri(x-\-api) -U^ ll (x)r\(x-api) where U-^(x) = llj l (x- fi) . (2.4) 
We have the correspondence 

15 3 7 

Gi <--,-,••■ and G 2 ^ (2.5) 

between the lattice and the D'^ continuum quantum numbers. For r > some 0' h /0(3)' 
irreps will split up into two or more D'^JD'^ irreps. For instance, the operators corresponding 
to H split up, depending on the direction, relative to the interquark axis: for the axis pointing 
into the 3-direction, we call the operator 71 Ai — 7 2 A 2 "parallel" (||) and the other combinations 
"perpendicular" (_L), where _L corresponds to the state of lower angular momentum (G\/\) and || 
to the higher angular momentum {Gj/\), relative to the 3-axis. 

3. Results 
3.1 Spectrum 

The spectrum of our observables at r = is shown in figure 2(a). The Regge trajectories 
are denoted by the grey lines and although they are not examined further they help to group the 
operators. The spectrum extends over a range of 2 GeV and it gets even richer when one goes to 
r > where the symmetry group 0' h breaks down into D' 4h . 

The states created by the operators P and P75 split up into two energy levels depending on 
whether or not the angular momentum contains a projection onto the QQ axis. The spectrum at 
r = 5a is shown in figure 2(b) and the grey lines in the plot show degeneracies of some of the 
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1.4 
1.3 
1.2 
1.1 
1 
0.9 
0.8 
0.7 
0.6 
0.5 



Masses of the operators at r = 
GS i- 




[GeV] 


[a] 


3.25 


1.8 


3 


1.7 


2.75 


1.6 


2.5 


1.5 




1.4 


2.25 


1.3 


2 


1.2 


1.75 


1.1 




1 
(1.!) 


1.5 


1.25 


0.8 



Masses of the operators at r = 5a = 0.42 fin 



Dir. 



D 



GS i — h- 

FE i 

SE jK- 



(a) r = 



P75|| -P75± P\\ 

(b) r = 5a 



P± 



2s 
75 



4.25 
4 

3.75 
3.5 
3.25 
3 

2.75 

2.5 

2.25 



Figure 2: The different operators at r = and r = 5a, with the lowest possible continuum spin assignments. 



levels. For example, the excitation of 1 and the ground state of D75 are degenerate. Both are in 

the G\ u representation of D' Ah , corresponding to the continuum D'^ quantum numbers | , | , 

The latter operator corresponds to the 0(3) quantum numbers | at r = 0. Hence we assign a 
continuum | spin to it and to the first radial excitation of the 1 operator. Unfortunately the two 
ground states of i , i.e. 1 and P _L seem to be different and it looks like \ ~ 1 is much closer to 

— ■ / / —II 

I ~ P || . We hope to resolve this issue in the near future. 




Figure 3: Potential for \ g ~ 750. Figure 4: Lowest lying potential for each operator. 



3.2 Potentials 

In figure 3 we plot the results of 5 ~ 75© as red, green and blue symbols. The factorization 
model expectation for the ground state m.Q q + \ Vqq is represented by the orange triangles. For 
r > 3a ?s 0.25 fm the data lie systematically above the expectation. The increase in energy with 
increasing r originates from the internal energy of the QQ diquark which, in the presence of the 
light quark, rises more steeply then expected. In figure 4 we see that this slope depends on the 
D' 4h quantum numbers, such that this failure of the factorization ansatz cannot be attributed to a 
different functional form of the QQ potential alone. 
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The pink and light blue lines are interpolations of the ground state points, shifted by the first 
two static-light energy splittings. The pink line describes the first excitation of the baryonic po- 
tential very well. This suggests that this excitation is due to the light quark, with very little effect 
from the gluonic flux configuration. The second excited state data lie somewhat above the light 
blue line. The most likely reason for this deviation is our inability to reliably isolate this excitation 
in our three-dimensional variational basis, so that we somewhat overestimate the masses. 

At large r one of the static quarks might form a diquark with the light quark. In this case we 
should see the transverse modes of the string connecting this Qq diquark with the remaining Q in 
the excitation spectrum. The Nambu-Goto string potential [13] suggests the following functional 
form: 

E ^ = GGsr \J l+ { 2n - d -w-)^- (3 - 1} 

In our case d = 4 and Cos is the effective string tension of the ground state determined by fitting to 

V QQq( r ) = Off + a eff r ~ ~ • ( 3 - 2 ) 

The red band is E2 —Eq + GS, which is the first candidate for an excitation because n = 1 would 
affect the spin of the system. This curve is above our first radial excitation and its shape is very 
different from the second excitation. The light quark excitation energies are smaller than those 
required for transverse excitations of the flux tube. The plots for the other observables look very 
similar. 

3.3 Wavefunctions 

We compute the RMS and the wavefunctions of the ground state and the first two excitations 
as outlined in [14]. The light quark RMS for the ground state wavef unction created by the operator 
75 is about 4.1a. Within errors it is independent of the geometry of the correlator. This distance is 
indeed reasonably close to r « 0.3 fm, where the factorization ansatz starts to fail. The shapes of the 
Coulomb gauge wavefunctions agree with the findings of [14] although this reference considered 
mesons while we study the light quark distribution inside a <2<2<?-baryon. 




Figure 5: Wavefunction of groundstate and first two excitations of 75© for r = 4. 



4. Conclusion and Outlook 

The factorization model breaks down for separations greater than 0.3 fm. The nature of the 
first excitation of a static-static-light baryon does not change with an increase of the separation in 
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the two static quarks. The second excitation could not be resolved so well and we do not venture to 
explain its nature. The splittings induced by breaking the symmetry group Oh into D4/, are reflected 
in our correlators. 

It is a puzzle why the states of \ u coming from the operators 1 and from Pj_ do not agree but instead 
the 1 seems to be degenerate with the | from Pn. At the moment we are conducting a study on a 
bigger lattice (24 3 x 48) at smaller quark masses (K = 0. 1362) and we hope to clarify this point. 
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